Role of the singular factors in the standard fits for initial parton densities 
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, Total resummation of double- and single- logarithms of x contributing to the spin-dependent struc- 

f*"^ ■ ture function gi ensures its steep rise at small x. In the asymptotic limit x — > 0, the resummation 

leads to the Regge behavior of gi and allows to calculate the non-singlet and singlet intercepts of g%. 
. DGLAP lacks such a resummation but suggests special phenomenological fits for the initial parton 

densities such that the singular factors x~ a in the fits mimic the resummation and also provide 
^ ' gi with the steep (power-like) rise at the small- a; region. Accounting for the total resummaton of 

h^- , logarithms of x allows to drop the singular factors in the fits and leads to a remarkable simplification 

^ . of the fits. 
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PACS numbers: 12.38.Cy 



J> 1 I. INTRODUCTION 

m ' 

The standard theoretical instrument for investigating the DIS structure function gi(x, Q 2 ) is DGLAP pj. In this 
approach, gDGLAP j g a conv olution of the coefficient functions Cdglap and evolved parton distributions which are 
also expressed as a convolution of the splitting functions Pdglap and initial parton densities. The latter are found 
. with fitting experimental data at x close 1 and Q 2 ~ 1 GeV 2 . There is an obvious asymmetry in treating Q 2 - 
and x- contributions in DGLAP. Indeed, the leading Q 2 - contributions, ln(Q 2 ), are accounted to all orders in a s 
whereas Cdglap(x) and Pdglap fa) are known in first two orders of the perturbative QCD. The reason for such 
O |' asymmetry is the fact that originally DGLAP was constructed for operating at large Q 2 and x not so far from 1. 
In this region ln fc Q 2 were large whereas x- contributions from higher loops were small and could be neglected. On 
the contrary, in the small- a; region the situation looks opposite: logarithms of x, namely double logarithmic (DL), 
i.e. the terms (a s ln 2 (l/x)) k 1 and single logarithms (SL), the terms (a s ln(l/a:)) fe ,with k = 1,2,.., are becoming 
quite sizable and should be accounted to all orders in a s . When the total resummation of DL terms was done[J 
it led to new expressions, g® L , for g\. In particular, Refs. Q demonstrated that the small- a; asymptotics of <?f L 
was of the Regge (power-like) form and was much greater than the well-known small- a; asymptotics of g p>GLAP (see 
;_i " Eq. at x < 0.01. However, this result was strongly criticized (see e.g. Refs. 0) where it was shown that using the 
standard NLO DGLAP formulae together with appropriate fits for initial parton densities led to the opposite result: 
gi L -C g p>GLAP at small x. After that the total resummation of logarithms of x for g\ was claimed irrelevant for 
available range of x. 

Strictly speaking, comparison of results of Refs. and could not be done in a straightforward way because 
DGLAP -evolution equations have always used the running a s , with the parametrization 

„ dglap „ ir&\ m 
a s = a s (Q ), (1) 

whereas Refs. Q operated with a s fixed at an unknown scale. Trying to set a scale for a s when the approach of 
Refs. is used, Refs. Q| suggested that the argument of a s should be Q 2 like in DGLAP. The parametrization 

a DGLAP = a 

s{Q 2 ) appears when the argument of a s in each of the ladder rungs of the involved Feynman graphs is 
kj_ , with k being the momentum of the upper parton (a quark or a gluon) of the rung. A deeper investigation of this 
matter Q led us to the conclusion that the DGLAP- parametrization of Eq. can be a good approximation at x not 
far from I only. Instead of it, a new parametrization was suggested where the argument of a s in every ladder rungs is 
the virtuality of the horizontal gluon (see Ref. for detail) . This parametrization is universally good for both small 
x and large x. It converges to the DGLAP- parametrization at large x but differs from it at small x. Using this new 
parametrization allowed us to obtain in Refs. |7j expressions for gi accounting for all-order resuummations of DL and 
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SL terms, including the running a s effects . In the first place, it was used to obtain numerical values of the intercepts 
of the singlet and non-singlet gi . It is worth to mention that these results were immediately confirmed 8] by several 
independent groups who fitted available HERMES data and extrapolated them into the asymptotic region x — > 0. 

Nevertheless, it is known that, despite DGLAP lacks the total resummation of In x, it successfully operates at x <C 1. 
As a result, the common opinion was formed that not only the total resummation of DL contributions in Refs. Q 
but also the much more accurate calculations performed in Refs. Q should be out of use at available x and might be 
of some importance in a distant future at extremely small x. In Ref. we argued against such a point of view and 
explained why DGLAP can be so successful at small x: in order to be able to describe the available experimental data, 
DGLAP uses the singular fits (see for example Refs. for the initial parton densities. Singular factors (i.e. 

the factors which — > oo when x — > ) in the fits mimic the total resummaton of Refs. Q. Using the results of Ref Q 
allows to simplify the quite complicated structure of the standard DGLAP fits down to normalization constants at 
small x. 



II. DIFFERENCE BETWEEN DGLAP AND OUR APPROACH 

In DGLAP, gi is expressed through convolutions of the coefficient functions and evolved parton distributions. As 
convolutions look simpler in terms of integral transforms, it is convenient to represent g± in the form of the Merlin 
integral. For example, the non-singlet component of gi can be represented as follows: 



9idglap( x i Q Z ) 



K/2) 



— {l/xyC D GLAp{u)5q{u) exp 
Ztir 



dk 2 

-prlDGLAp{u, a s {kj_)) 

Kj_ 



(2) 



with Cdglap{u) being the non-singlet coefficient functions, 7dglap(w, a s ) the non-singlet anomalous dimensions 
and Sq(uj) the initial non-singlet quark densities in the Mellin (momentum) space. The expression for the singlet g\ 
is similar, though more involved. Both "/dglap and Cdglap are known in first two orders of the perturbative QCD. 
Technically, it is simpler to calculate them at integer values of u> = n. In this case, the integrand of Eq. is called 
the n-th momentum of g^ s . When the moments for different n are known, g NS at arbitrary values of u> is obtained 
with interpolation of the moments. Expressions for the initial quark densities are defined from phenomenological 
consideration, with fitting experimental data at x ~ 1. Eq. (J2J) shows that jdglap govern the Q 2 - evolution whereas 
Cdglap evolve Sq(oj) in the x-space from x ~ 1 into the smalkc region. When, at the cc-space, the initial parton 
distributions Sq(x) are regular in x, i.e. do not — > oo when x — > 0, the small- a; asymptotics of g\ dglap is given by 
the well-known expression: 



9l DGLAP' 9l DGLAP ~ ex P 



Hl/x)ln(HQ 2 /^ 2 )/Hv 2 /A 2 QCD ) 



(3) 



On the contrary, when the total resummation of the double- logarithms (DL) and single- logarithms of x is done 5], 
the Mellin representation for g^ s is 



g? s (x, Q 2 ) = (e 2 /2) \ ^L{l/ X yC NS {u)8q{u>) exp {H NS {u>) \n(Q 2 /^ 2 )) , 



(4) 



with new coefficient functions Cats, 



Cjvs(w) 



and anomalous dimensions H^s, 



where 



w - H NS (u>) 



H NS = (1/2) |w - y/u; 2 - B(lo) 
B(lu) = (4ttCf(1 + lu/2)A(uj) + D(lu))/(2it 2 ) 



(5) 



(0) 



(7) 



The parametrization of see Ref. |5| was used later in Refs. ;@j for studying the small-x contribution to the Bjorken sum rule. 
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D(ui) and A(u) in Eq. JJJ are expressed in terms of p = m(l/x), r\ = m(/i 2 /Ag C£) ), = (33 — 2n/)/127r and the color 
factors C F = 4/3, N = 3: 



2C F 



d j oe-" p In ( 



P + V 



{p + rif 



1 



p + t] 



(8) 



'7 



?/ 2 + 7T 2 



Hs and Cats account for DL and SL contributions to all orders in a s . When x — » 0, 



ft™ 



(JfC?)*""' 2 , gf 



(9) 



(10) 



where the non-singlet and singlet intercepts are Ajvs = 0.42, A5 = 0.86. The x- behavior of Eq. (|10fl is much steeper 
than the one of Eq. ©. Obviously, the total resummation of logarithms of x leads to the faster growth of g\ when x 
decreasing compared to the one predicted by DGLAP, providing the input Sq in Eq. @ is a regular function of to at 

III. STRUCTURE OF THE STANDARD DGLAP FITS 

Although there are different fits for Sq(x) in literature, all available fits include both regular and singular factors 
when x — * (see e.g. Refs. 0,^1 for detail). For example, one of the fits from Ref. i s 



Sq(x) = Nr}X~ a (1 - xf {I + 1X °) 



(11) 



with N, 77 being a normalization, a — 0.576, (3 = 2.67, 7 = 34.36 and 6 = 0.75. In the 10 -space Eq. is a sum of 
pole contributions: 



Sq(uj) — Ni] (uj — a) 1 + ^ m k (uj + X k ) 



(12) 



with Xk > 0, so that the first term in Eq. I)12|l corresponds to the singular factor x~ a of Eq. 111(1 . When the fit Eq. 1|11|) 
is substituted in Eq. J5J), the singular factor x~ a affects the small -x behavior of 171 and changes its asymptotics Eq. © 
for gi for the Regge asymptotics. Indeed, the small- x asymptotics is governed by the leading singularity u — a, so 



gi DGLAP 



C(a)(l/x) a [(HQ 2 /A 2 ))/(H^/A 2 )) 



■y(a) 



(13) 



Obviously, the actual DGLAP asympotics of Eq. i|13|) is of the Regge type, it differs a lot from the conventional 
DGLAP asympotics of Eq. © and looks similar to our asymptotics given by Eq. (|10|) : incorporating the singular 
factors into DGLAP fits ensures the steep rise of g^ GLAP at small x and thereby leads to the success of DGLAP at 
small x. Ref. demonstrates that without the singular factor x~~ a in the fit of Eq. (|llf> . DGLAP would not be able 
to operate successfully at x < 0.05. In other words, the singular factors in DGLAP fits mimic the total resummation 
of logarithms of x of Eqs. (|4I10() . Although both (|13fl and (|10fl predict the Regge asymptotics for gi, there is a certain 
difference between them: Eq. (|13f) predicts that the intercept of g^ s should be a = 0.57. As a is greater than the 
non-singlet intercept A^s = 0.42, the non-singlet g DGLAP grows, when x — > 0, faster than our predictions. However, 
such a rise is too steep. It contradicts the results obtained in Refs. and confirmed in Refs. Usually, the DGLAP 
equations for the non-singlets are written in the x-space as convolutions of splitting functions P qq with evolved parton 
distributions Aq and the latter are written as another convolution: 



Aq(x) = C q (x,y) ® Sq(y), 



(14) 



with C q being the coefficient function. Written in this way, Aq is sometimes believed to be less singular than Sq 
because of the evolution. However applying the Mellin transform to Eq. l|14fl immediately disproves it. 
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IV. CONCLUSION 

Comparison of Eqs. 10 and l|13|) shows explicitly that the singular factor x~ a in the Eq. Ullfl for the initial quark 
density converts the exponential DGLAP-asympotics into the Regge one. On the other hand, comparison of Eqs. i|l(J|) 
and (|13|l demonstrates that the singular factors in the DGLAP fits mimic the total resummation of logarithms of 
x. These factors can be dropped when the total resummation of logarithms of x performed in Ref. Q is taken into 
account. The remaining, regular ir-terms of the DGLAP fits (the terms in squared brackets in Eq. (|llll ) can obviously 
be simplified or even dropped at small x so that the rather complicated DGLAP fits can be replaced by constants. 
It immediately leads to an interesting conclusion: the DGLAP fits for Sq have been commonly believed to represent 
non-perturbative QCD effects but they actually mimic the contributions of the perturbative QCD, so the whole impact 
of the non-perturbative QCD on g\ at small x is not large and can be approximated by normalization constants. 
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